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ABSTRACT 


The  heating  of  ions  in  a  magnetized  plasma  by  the  use  of 
the  second  order  fields  generated  by  the  nonlinear  mixing  of 
two  high  frequency  waves  (  high  frequency  means  U>  y  U3f>_  ) . 

The  kinetic  equations  describing  the  mixing  and  heating  process 
are  solved  using  the  method  of  orbit  integrations.  The  energy 
is  absorbed  from  the  second  order  fields  via  cyclotron  damping. 
The  sensitivity  of  the  absorbed  power  to  fluctuations  in 
frequency,  density,  static  magnetic  field  and  mixing  angle 
is  analysed. 
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CHAPTER  1  INTRODUCTION 

Controlled  thermonuclear  fusion  has  been  on  the  minds  of 
engineers  and  physicists  since  the  explosion  of  the  first 
hydrogen  bomb.  A  large  amount  of  time  and  money  has  been  spent 
in  the  last  two  decades  to  achieve  this  end.  Although  consider 
able  progress  has  been  made  towards  achieving  this  goal,  the 
basic  problems  still  to  be  surmounted  are  essentially  the  same 
as  they  were  at  the  onset  of  research.  The  two  basic  problems 
associated  with  fusion  are  confinement  and  heating. 

Since  this  thesis  deals  only  with  plasma  heating,  a  dis¬ 
cussion  of  the  state  of  the  research  in  plasma  confinement 
will  not  be  given.  (For  a  detailed  discussion  of  plasma  con¬ 
finement  techniques,  see  references  1-4.) 

To  achieve  positive  energy  yields  from  a  fusion  reaction, 

5  8 

it  is  estimated  that  ion  temperatures  of  IO  °Kor  better 

will  be  needed.  A  first-order  calculation  of  electron-elec- 

6 

tron  and  electron-ion  collision  frequencies  done  by  Spitzer 
gives 


(1.1) 


$  A  =  (kT 

\  e.  • 
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and 


(1.2) 


Equations  1.1  and  1.2  show  that  the  collision  frequencies  in 


a  nlasma  are  inversely  proportional  to  T  L 


This  in  turn 


implies  that  any  heating  method  which  depends  on  electron- 
electron  collisions  for  energy  absorption  will  become  more 
and  more  inefficient  as  the  temperature  of  the  electrons  in¬ 
creases.  Equation  1.2  indicates  that  ion  heating  via  the 
equipartition  of  energy  from  the  electrons  is  a  poor  mechanism 
at  high  temperatures.  From  the  conservation  of  momentum, 

.  For  the  case  where  W*.  (required 

for  charge  neutrality  in  a  two  component  plasma). 


Equation  1.3  implies  that  if  the  ions  are  selectively  heated 
(ion-electron  equipartition  times  are  Yt\*t  |we  larger  than 
electron-ion  equipartition  times)  a  significant  amount  of 
electron-ion  temperature  decoupling  can  occur.  This  means 
that  the  ions  can  in  principle  be  maintained  at  a  higher  tem¬ 
perature  than  the  electrons.  This  propertv  of  a  plasma  is 
very  important  since  in  order  to  achieve  fusion  one  should 
"ideally"  have  ions  at  lO  °Kand  cold  electrons.  This  is 
because  the  electrons,  due  to  their  high  mobility,  are  responsible 
for  the  majority  of  the  losses  from  a  hot  plasma.  The  two  most 
important  mechanisms  for  losses  via  electrons  are  cyclotron 
and  bremss trahlung  radiation.  It  is  clear  from  this  argument 


■*, 


, 
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that  the  most  efficient  heating  methods  to  be  emoloved  for 
fusion  will  be  those  which  selectively  heat  the  ions  through 
a  non-collisional  mechanism.  In  order  to  achieve  this  end  a 
considerable  amount  of  research  in  the  field  of  nlasma  heating 
has  been  directed  at  the  problem  of  ion  energy  absorption 
via  cyclotron  damping. 

Heating  Techniques  Proposed  to  Date 

In  this  section  a  brief  review  of  the  various  methods 
used  to  date  for  plasma  heating  is  given. 

1.  Ohmic  Heating 

Plasma  Ohmic  heating  using  high  intensity  laser  radiation 

has  recently  been  proposed  as  a  method  by  which  a  plasma  can 

be  brought  to  thermonuclear  temperatures.  Ohmic  heating  has 

7 

been  used  in  fusion  machines  for  preheating  the  plasma.  Ohmic 

heating  arises  strictly  from  the  resistivity  of  the  plasma 

to  the  electric  field.  Since  the  collision  frequency  varies 
-3/2_ 

as  T  '  ,  this  method  of  heating  becomes  more  and  more 
inefficient  as  the  electron  temperature  increases.  Another 
disadvantage  of  this  method  is  that  most  of  the  power  is  ab¬ 
sorbed  by  electrons.  At  high  temperatures,  the  electron-ion 
equioartition  time  becomes  large,  making  the  transfer  of  energy 
between  the  electrons  and  ions  very  inefficient. 

2.  Magnetic  Pumoing 

8,9,10 

In  this  scheme,  the  magnetic  field  confining 
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the  plasma  is  modulated  in  time.  Heating  in  the  perpendicular 

direction  occurs  because  of  the  adiabaticity  of  the  magnetic 

moment  and  the  absorbed  energy  is  transferred  via  collisions 

8,9 

to  the  parallel  direction.  This  technique  involves  the  ion- 
ion  relaxation  time  and  is  therefore  very  inefficient  at  high 
temperatures.  It  also  suffers  from  the  fact  that  at  high 
frequencies  (where  the  heating  would  presumably  be  more  ef¬ 
ficient)  only  a  very  thin  surface  layer  of  the  plasma  is  com¬ 
pressed  and  only  the  surface  layer  is  being  heated. 

3.  Transit-Time  Heating 

In  transit-time  heating,  the  magnetic  field  is  modulated 

8 

in  time  and  space.  In  this  method,  the  longitudinal  adiabatic 
invariant  is  responsible  for  energy  absorption  and  the  ran¬ 
domization  occurs  via  large  scale  scattering.  A  serious  draw¬ 
back  to  this  method  is  the  splitting  up  of  the  original 
Maxwellian  distribution  into  a  fast  and  slow  component  leading 
to  instabilities  in  velocity  space.  Numerous  theoretical 
and  experimental  investigations  carried  out  in  the  last  few 
years  on  transit-time  heating  have  shown  this  scheme  to  be  the 

9 

strongest  contender  for  heating  plasmas  in  toroidal  geometries. 

A  more  detailed  discussion  of  magnetic  pumping  can  be  found 

11 


in  James  and  Capjack. 


- 
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4.  Rapid  Compression 

A  preheated  plasma  column  may  be  heated  through  a  radial 

compression  obtained  by  pulsing  up  the  strength  of  the  magnetic 

field.  The  velocity  imparted  to  the  particles  using  rapid 

compression  is  species  independent.  This  species  independence 

implies  that  the  ions  are  gaining  most  of  the  energy  from  the 

compression.  An  experiment  utilizing  this  type  of  heating 

scheme  is  the  2X  mirror  experiment  at  the  Lawrence  Radiation 
12 

Laboratory . 

5.  Ion  Cyclotron  Resonance  Heating  (I.C.R.H.) 

Plasma  heating  by  the  use  of  frequencies  at  or  near  the 

ion  cyclotron  frequency  has  been  under  intense  theoretical 

13-17 

and  experimental  investigation. 

a.  Reverse  Turn  Induction  Coils 

The  coil  is  often  referred  to  as  a  Stix  coil  and  is  used 

to  generate  cyclotron  waves  in  a  plasma  loaded  cylindrical 
18 

wave  guide.  The  Stix  coil  configuration  has  been  used  in 

the  B-65  stellerator  at  Princeton  and  has  also  been  used  in 

19 

the  model  C-stellerator .  The  coil  is  built  so  that  current 
passing  through  a  section  is  180°  out  of  phase  with  the  current 
in  the  adjacent  sections.  This  allows  the  electrons  to  flow 
along  the  magnetic  field  lines  to  cancel  the  radial  oscillating 
space  charge  electric  field  that  would  otherwise  be  set-up, 
resulting  in  the  plasma  shielding  itself  from  the  applied 


field. 
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b.  Coaxial  Electrodes 

R.F.  power  is  applied  to  coaxial  electrodes  of  a  mirror 
device  and  a  torsional  Alfven  wave  is  generated  in  the  plasma. 
Since  the  electric  field  of  the  wave  in  the  plasma  is  in  the 
radial  direction,  an  efficient  coupling  of  energy  is  possible. 
The  theoretical  coupling  efficiency  can  be  estimated  by  treating 
the  device  as  a  coaxial  wave  guide  filled  with  a  material  of 
dielectric  constant  \  The  waves  are  generated 

in  a  region  where  « 5  Absorption  of  the  wave  energy 

is  then  accomplished  through  the  use  of  a  magnetic  beach. 

20 

In  an  experiment  performed  by  Boley,  Wilcox,  et  al, 
the  transfer  of  energy  from  a  1MW,  8.3  MHz  oscillator  to  a 
plasma  with  a  density  of  6  x  10^  Cm-^  was  at  an  efficiency 
of  65%.  Better  than  90%  of  the  wave  energy  was  absorbed  by 
the  ions  at  the  magnetic  beach. 

21,22 

A  similar  experiment  was  performed  by  Shvets  et  al 
on  the  Vikhr  (Whirlwind)  device.  In  this  experiment  an  ion 
temperature  of  2.50  e.y.  was  obtained  in  a  plasma  with  a  density 
2  x  10^-4  cnf^. 

c.  Difference  Frequency  Ion  Heating 

Since  the  topic  of  this  thesis  is  plasma  heating  via  non¬ 
linear  mixing  of  waves  in  a  plasma,  a  more  detailed  discussion 
of  the  history  of  this  method  will  be  given. 


* 


A  large  amount  of  theoretical  and  experimental  work  has 

23-35 

been  done  on  non-linear  mixing  of  waves  in  plasmas. 

36 

In  1966,  James  and  Thompson  proposed  using  the  difference 
frequency  harmonic  which  is  generated  by  the  non-linear  mixing 
of  two  high  frequency  waves  for  plasma  heating.  (high  fre¬ 
quency  waves  are  defined  as  to  ensure  penetration 

of  the  incident  waves  into  the  bulk  of  the  plasma.)  In  order 
to  enhance  the  second-order  fields,  they  tuned  their  difference 
frequency  to  a  natural  resonance  in  the  plasma.  Using  a  plasma 
natural  resonance  made  the  absorbed  power  critically  dependent 
on  plasma  parameters  such  as  number  density  and  magnetic  field 
strength.  Stability  of  plasma  parameters  of  better  than  one 
part  in  10^  were  required  in  order  to  maintain  the  absorbed 
power  at  an  acceptable  level.  They  also  found  that  the  power 
absorbed  from  the  difference  frequency  harmonic  varied  as  \£,\ 
This  meant  that  for  the  t>ower  absorbed  would  be  very 

small  unless  the  electric  fields  were  very  large.  In  order  to 
lower  the  incident  frequencies  and  still  ensure  good  penetration 

of  the  incident  waves  in  the  plasma,  an  analysis  was  done  by 

37 

Capjack  and  James  using  incident  fields  in  the  whistler 
regime.  In  this  analysis,  the  difference  frequency  wave  was 
tuned  to  the  ion  cyclotron  resonance  and  a  full  kinetic 
analysis  was  done  to  calculate  the  power  absorbed  by  cyclotron 
damping.  This  analysis  also  proved  to  have  very  high  require¬ 
ments  on  the  static  fields  and  the  difference  frequency. 


■ 
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In  order  to  circumvent  these  problems,  a  kinetic  analysis 
of  cyclotron  damping  of  high  frequency  waves  is  undertaken  in 
this  thesis.  The  incident  frequencies  used  in  the  numerical 
calculations  are  chosen  such  that  a  CO^  laser  can  be  used 
to  experimentally  verify  the  results  of  this  analysis.  The 
difference  frequency  wave  is  set  equal  to  the  ion  cyclotron 
frequency.  The  power  absorbed  by  the  ions  and  electrons  at 
the  difference  frequency  is  computed.  Cyclotron  energy  ab¬ 
sorption  proves  to  be  relatively  insensitive  to  magnetic  field, 
difference  frequency  and  number  density  fluctuations.  This 
method  of  energy  absorption  shows  an  improvement  as  the  number 
density  of  the  plasma  increases.  These  characteristics  make 
this  method  a  good  candidate  for  heating  to  thermonuclear 
temperatures . 
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CHAPTER  2  NONLINEAR  INTERACTION  OF  WAVES  IN  A  PLASMA 


In  this  chapter  the  collisionless  Boltzmann  equation 
will  be  solved  to  second  order  to  describe  the  nonlinear  mixing 
of  two  high  frequency  waves  in  a  plasma.  To  solve  the  Boltzmann 

38 

equation  the  method  of  orbit  integration,  first  used  by  Drummond 

to  calculate  first  order  fields  and  currents  and  latter  extended 

37 

by  Capjack  and  James  to  include  second  order  effects,  will 
be  used. 


2.1  Kinetic  Equation 

Consider  a  particle  trajectory  given  by 

$  = 

Then  the  total  rate  of  change  of  the  j  '  th  distribution  along 
this  trajectory  is  given  in  the  Lagrangian  coordinates  along 
the  trajectory  by 

=  iSi  *  *■  4v 

•at  At  3*  At 


At 


(2.1) 


Setting  Vs  ,  using  the  Lorentz  force  equation  to  replace 


,  and  taking  the  limit 


^  by  v 

At  WJ  c 

of  no  collisions  (  ~~  O  )  equation  2.1  becomes  the  Boltzmann- 

Vlasov  equation. 

w’ 


(2.2) 
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The  Vlasov  equation  coupled  with  Maxwell’s  equations  for  the  E.M. 
field  V  X  E  -  -X^£> 

VxB=  V*E.=  <RP  (2.3) 

C3t  1 

and  with  the  definition  of  charge  and  current 

/>-  §e^£  SE&A 

form  a  closed  set  of  equations  which  can,  in  principle,  be  solved 
f°r  . 

£note:  The  validity  of  using  the  Boltzmann-Vlasov  equation  is 
discussed  in  standard  books  of  Plasma  Physics  (i.e.  Gartenhaus- 
Elements  of  Plasma  Physics-chapter  4)  and  will  not  be  discussed 
here.  It  will  also  be  implicitely  assumed  that  a  second  order 
perturbation  expansion  of  T.  provides  an  "adequate"  description 

b 

of  wave  mixing  in  a  hot  magneto-plasma!} 

2.2  Solution  of  the  Kinetic  Equation 

The  solution  to  the  set  of  equations  2.2  and  2.3  is  taken  to 


be  an  expansion  to  second-order  of  the  form: 


* 


4©  4m  -h> 
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‘-(fc^S  -<aO| 
e  + 


iL  tSfX  - 
e 


it 

5J 


U)C) 

e  +- 


complex  conjugate  terms 


(2.4) 


where  K-  K,-Kt  and  VO  -  VO^  -VOj, 

The  £  signs  refer  to  the  sign  of  the  charge  of  the  two  component 

plasma  under  consideration.  The$"  are  the  zeroeth, 

TO  >  Tl,7.,  TS 

first  and  second  order  terms  to  the  distribution  function.  To  make  the 


mathematical  analysis  of  this  problem  tractable  it  will  be 
assumed  that 


Bo-  B.2 

4  (2.5) 

E,  =  4  v 

Ejl  a  ^5t 
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2.2.1  Solution  of  the  Zeroeth-Order  Equation 

The  zeroeth-order  trajectory  is  defined  by  the  equation 
of  motion 


*  %  z-ee \/xb„ 

tit  &  ^  " 

The  rate  of  change  of  X  along  the  zeroeth  order  trajectory 
is  given  by 

(ii\  =  v 

M+  'o  -at  ^  Ytl-t 

The  rate  of  change  of  the  zeroeth-order  solution  of  the  B.V. 


(2.6) 


(2.7) 


equation  along  the  zeroeth-order  trajectory  is 

(41A  *  -v  N  •  Zee  O  (2. 7A) 

' dv  K  -at  a£  mr  "  ^ 

from  the  B.V.  assumption.  The  time  and  space  independent 

solution  of  2.7A  is  the  form  tM.H'O*')  where 
refer  to  directions  with  respect  to  the  static  magnetic  field. 

In  this  analysis  it  will  be  assumed  that  the  zeroeth-order 


distribution  is  Maxwellian. 

(a!-  f 

°  VfcTTKT*) 


l.  e. 


e  z*rt 


(2.8) 


2.2.2  First-Order  Solution 

The  first  order  solution  to  the  Boltzmann  equation  is 
obtained  by  expanding  to  first  order: 

r  =  £  *  £; 

Using  this  expansion  in  2.2  and  retaining  terms  of  first  order 
yields: 


^  ti  4-  gee  ( 6^  •  Ze€.(£t^^x.sh-^jf ) 

mt  *•  3  V  tnt 


. 
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If  the  total  derivative  of  is  taken  along  the  zeroeth- 
order  trajectory  and  the  result  substituted  in  2.9,  then  2.9 


becomes : 


Rf*)  - 


At  \  XYlt 

This  equation  can  be  integrated  along  the  zeroeth-order 
trajectory  to  yield 

I,1  =  -  \  (i,1  +•  V  | ^ 

va  v  c- 

[Note;  Rigorously  speaking,  the  integration  should  have  been 


(2.10) 


written: 


^  =  -lee  +  11'  +  |1(XX10 

TO  A  ^  1 


However,  if  attention  is  focused  only  on  solutions  that  are 
growing  in  time;  then  4  becomes  exponentially 

small  as  4*  — —  OO 


and  may  therefore  be  neglected.] 


2.3  Second-Order  Distribution  Function 

A  second-order  solution  to  the  B.V.  equation  is  obtained 
by  using  the  expansions  2.4  in  equation  2.2.  .  /  u  +  .  %i.\ 


.  r  t  Kt  itK.r-w.O  t  ‘•(Kt-r-vuiV>  t 


+  c 


.C.1+ 

4  1 7 


Bo+  §ie 


LCK,r-vj,n  iCwr-vo^ 

+•  ; 


fc  c  c.J- i_fr- 

u.e  +  £*&  4t*e  ’  J  av  \y» 


* 


*  ^  UXz-Sr'^  ftt 


ti'  e 


i 


+  t3e 


4-  c.c* 


■] 
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Using  harmonic  balance  to  equate  terms  with  common  exponential 

C (K -X 

^  where  K55 


and  CO  a  vxi  ^ 


yields : 

2^5T  +  v.3^s  +  Zee(vxft  = 


•at 


ar 


xnt  c  3^ 


-  iee 


c 


(2.12) 


(e,  i. x*_0 * 1 

c  n  J 


3* 


Since  along  the  zeroeth-order  trajectory 


£  dJ.UIIg  Lilt:  4CLUeLll-ULUCl  LldJCLLUiy  ^ 

MS^N  =  3|?  +  ^-^i  +  zeg  Cv^oV^' 

1  Xt  )n  at  3t  Wt  *  at 


(2.13) 


then 


(— 3  \  =  -zee  {(e3  +  vxgjV^JcT  +  ( +  Vxft^  •  +■  l 

MV  '  YYW 


By  integrating  equation  2.14  along  the  zeroeth-order  trajectory 


^  f  ( c3' t  sV&ValT  At' 

mt  J  -=  3v' 

-oO  ^ 

s  r  c  £,* + x,£§t'v  At 

*.i  c  ** 

ffc  «.*• 

-Jco  c  av' 


-  ~£ee 
rat 


-  g.ee 


(2.15) 


* 
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t' 


can  then  be  written  as 


where 


V  ■  K  *  &  * 

•i*  »  -uee  (  t  At 

Tl*  m*  J  c.  av* 


(2.16) 


(2.17) 


fu  a  "  ^  \  ( c«'  *  ^ '  =L?a  <&' 

_Jte  c  ay 


f  (u;* 

Ta.  rnt  e. 


(2.18) 


(2.19) 


Note:  Specific  reference  to  the  species  of  the  particle  will 

be  dropped  at  this  point  since  necessary  signs  and  factors 
are  conserved  by  terms  like  <1^  , 


2.3  Evaluation  of  X 

Til  f 

To  evaluate  an  expression  for  is  required.  This 

is  obtained  by  using  the  Fourier  analysed  Maxwell's  equation 

Sfil  ^1  (2.20) 

to  replace  by  %  i  in  2.10  #to  get 

4's  -uec  ( i  +  (2.2i) 

*  «t  -«  7x7  ui,  ax' 

For  Maxwellian,  2.21  reduces  to  ****  ** 


i: 

^  mm  (A 


Wl 


(2.22) 
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The  integrand  of  2.22  must  be  evaluated  along  the  zeroeth- 
order  trajectory  defined  by  r  =  X  (V* 

.  The  solution  to 

2.6  which  reaches  Va-Y(  at  is 

yx"  =  vx  -  tVvj'SiwnCt'-f) 

S'  =  evx'  +v^  Co^nCt'-t"') 

v2‘ 


(2.23) 


Th 


e  integration  of  2.23  yields 


t  tW 


x"=  -V  S^n.(t'-t,N)  +  ev^1 1 1 -  Coi &CV-t“)}  +-x' 

■a  o. 

u>"  =  -£Vxti-CcA Si.Cjf  -t"li  -  Hi  Sk«Sl tt  -tuj  +  vj  i 

£1  Si  J  (2.24) 

1."  =  -Vs-'Ct’-t'O  +  z/ 


By  using  2.23  and  2.24  to  replace  all  doubled  primed  quantities 


by  primed  quantities 


t:- 


L  (a,vx'  +  V>,vJ  +  4- 


Lfa-r'-uX)  l  Vu,,x  T  D'") 
gee  e  Hv4  \  e  {(y*  CosaM  - 
1  va 


+  I  At 


(2.25) 


Where 


ai-  -Kx,^ax 


?  =  t'-t" 


(2.26) 


.  - 
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b»  -  K ac.t  L  V  -  Cos-Sl/tl 

vSL 


(2.27) 


C1- 

(2.28) 

i 

By  using  equations  2.20  to  express  the  magnetic  field  in  terms 
of  the  electric  field  in  equation  2.19 

^  J  6  I  (l- V  A.t 
J|  m  -oo  V.  BVi  J 

and  using  2.25  for 


(2.29) 


■I  — ^ 

^31  -* — ; 


e  |  (Vx 

-«*  o  V. 


xCo53101  -  eV^  SmSt/fc 


1  ~  ^jgVz  X  lQ>i  ~  •**  Vx.  CcC,  -  vtW-*  (2.30) 

te,*  »oOx  ICT  Uiz  VCV  J 

+  (l- KrtV^CosSt-t  v  b,t  l 

£|X  ^2.  J 

The  integration  with  respect  to  "t  must  be  calculated  along  the 
zeroeth-order  tralectorv  defined  by  X.=jftO  .  The  solution  which 
reaches  V*=.  X '  V'Uy^at  t*  =b  is  obtained  from  2.23  and  2.24  bv 
replacing  "t  by  ”t  and  "tby  t  .  Employing  these  new  equations 
and  substituting  for  "fc—  t  yields  after  considerable 


algebra: 
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‘51  - 


r*°  r’°  <■  (  ■s'Xsc .  +  +  (u\  -KijVj')  ■* 

J  J  efA»+  &\*V*  +  6,,^ 

D  O  L 


e"  g,>  e  * 

lm\<T 


+  (ui-  K,V^,u' ) 

+  6,4  Vo*  +  01,  V,’  -V  4n,V*Vv|j  A'tck'X* 

The  &  S  can  be  expressed  in  the  more  useful  form  (using  ex¬ 
ponential  notation) 


(2.31) 


A 


u 


*,s- 


6,s  - 


ISlt"  -  i Sl't1  -Cslx"  +  i-Slt* 


°)n 


att=  e 


is-t'r 


-iSiZ-t":  \  -iQX'f  utiVN  \ 

V‘  +  e-  M*) 


ieSlt‘1  lesi'T  .  -I six' 

l  ^ise  ( 


icsi'r1 


(i 


H 

-c 

e 


c*  isi^N  \ 

-e  )v^) 


-Ua't‘7  isit"  -c»V'\  n 
+  e  Li%e  le  “« 


Bh  =  *n 


15 


(  ««=■(  ^'uw') . 


31A) 


fivu= 


/ -In**1  ^  Usir"  -ISLlf  ^  ^cS-T1* 

-cC-SlfM  -tCSTt 


') 
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where 


l(» 

2. 

***■  £\xS*. 

-x( 

-V 

2. 

101  £tx-£> 

x- 

-(»■ 

-  ^ 

A  SI 

( 

Z-Eix.  ^  ^ 

-3 

% 

11 

-lC 

W\  + 

4 

XT 

a  - 

-  \C* 

Swv^iC't.'V  ^ l>) 

w= 

'l.I-  Co^aUv't’i] 

XT 


(2.31B) 


2.4  Velocity  Moments  of  *3^ 

In  the  solution  to  this  problem,  the  quantities  of  interest 
are  the  zeroeth  and  first  moments  of  our  distribution  over 
velocity  space.  These  moments  are  of  the  form: 


oft  '-ft 


—  Oft  — o 


(2.32) 
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The  calculation  of  these  moments  involves  integrals  of  the 

form  i,  f  (cGN  -  WY\ 

'  -ZKT  ' 


G1  ^ 


Av 


which  can  be  evaluated  by  completing  the  square  to  yield 


-  a  icT 

G(o,aS  * 

"itn 

e 

_a?KT 

G(i,^  * 

Lator  e 

-a  KT 

GU.cft  = 

(  ut_  oTk'tA 

lYW 

e 

v  m  \nl  / 

i  a\<x(  3kT 

~  cf  k2tx 

m  ' 

Vf\' 

7.rt\ 


In  carrying  out  the  integration  over  Y*  and  one  obtains 

terms  in  the  exponent  of  the  form 

2.  TO 

which  can  be  evaluated  to  yield 


TO 


where 


Kx  XT 
si?-  to 


(2.33) 


(2.34) 


(2.35) 

(2.36) 


/  a»  *  i\  ■  *j  '  '  ! 
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The  quantity 


which  is  the  ratio  of  the  square  of  the  Larmor 


radius  to  the  perpendicular  wavelength  is  an  important  expansion 
parameter.  In  this  thesis,  consideration  will  be  given  mainly 
to  the  case  where  ^  is  large.  The  regime  of  large\  corresponds 
to  a  plasma  where  the  perpendicular  wavelength  of  the  difference 
frequency  is  much  smaller  than  the  electron  Larmor  radius.  With 
the  help  of  2.34,  one  can  evaluate  the  zeroeth  and  first  moments 


of 


e>0 


tnkT 


(2.37) 


t^ioAor  +  +  i3is f S - 

■WT  *\  Vn  1  v*1  ) 

J 


cC\VV'*^ 


<  ^  -  z*e*E  »C_e. 

i_H  -  v-  c»sn.(.,x+t''f) 

e 


ft  tw 


(2.38) 


f  B»  L\»\CT-  )  *  6,4 

[  "  »n* 


v  vw 


-  +  (VtWvcr/aw-  ♦  6(wt<^fsS-V5<V)\  ^kx' 

r*L  •  Vr\  \  Vrv  tn*  '  tn  W  h1  W 
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< 


K 


t^v 


&  ca\cT 
*n 


+  a^)  -  ^aWn'  *  o?jM 

vtf-  vir  ^  vw*  / 


+  &«  1^5 (a. «£!>')  ,  Aki^r(g-oagf)]i«it 


(2.39) 

I 


¥  =  (w-uy  t(w,-K„V*\X 


the  term  G.  can  be  expanded  by  using  the 


Bessel  function  equality 


S  e. 


Y\, 


(2.40) 


To  perform  the  integration  of  r^1  and  ,  one  must 

collect  common  exponential  terms  which  means  that  terms  like 
(XjId*  must  be  expanded. 


. 
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Ci 

e 

m 


C3*/tu  -cJl'X 

-  e 


cx^WY'j1  =  -<<*(  e 


c2A/t,‘ 


-IxS*.'*1* 
-  +  e 


) 


0?  h£EV  '  -uH  e  -  3e.  ■*•■!) 

'I'ft  I 


ed  ( 

(-)*  = 

<£( 

\\nn  * 

(KT\  * 

feol3 

^  Yt\  / 

(  z  - 


os\A 


-  e 


(2.41) 


13,,  „z  /  CS/tM  -uSLt«\  ,  /  -AzJlVl\ 

d  |  k^\  »  fed  V  Zo  -  l?(e,  -v  e.  )  V  lpi,e.  -v^  j 

i  .  r  f  /  iatM  -Xa*r*\  czs/e"  -czjiV1  \ 

b^VO^  -  <.6olH  l(-e  +e.  )  -  e  v  e  } 

,.  .,  H  ,(  iSW"  -Va**  f  ciAt”  \ 

oft  |UT\  =  -  e^H  -  4  4  e  +  e.  +  i^e.  +  <=,  ) 


El 


OStV 

“  e  -  e. 


) 


.  ~(  i  l5ttM  -lMw\  (  » 

r  UH SU  -e  )-AU 


4-e 


cWt*  -i-VSl^X 


") 


9  '  1  *\ 
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where 


ei  -  V^T 

za.  Vn 


Using  the  set  of  equations  2.41  for  the^'^V's  etc.,  2.40 
for  the  expansion  of  C*  ,  equations  2.37-2.39  can 


be  written  as 


c(!£<-uA^  ^  rr  l*\  +  ” 


<  1>A  *  « 

t*WT 


s  n 


^1 2. 


vSi/t1 


-I  si?*  t  r  ^  _  <siV  /-r-  — 


-Asi*' 


-e((  z(X^_l+€-X 


v\vv^ 


”  "^*V\-2L+€.  ^  -^V\.V2.+C  ^  ^ 


I 


(  2.  (  -  Xn-|_^  *  ^n-2.-e  “ 


^  +  c*v(. 


(wt  •■^YwOe*’**  +•  (rA.(  -IA+^e  ^  +  fe  (  (  tx„-  XAH-XA_^ 


le. 


-  Ut,-iw.xMVw-'v?((( 


nn 


UiV 


*  +  2.x e  +  C C  f^-0  X 


v\H 


v  ‘vZ'^-H*Oe  *  2,6  ^  £•("£*  -i-t 


^’A-'Z.-C  ^ 


U«si*' 


-  C  2’^^*A-VV€.  “"^n-H-6^  ~  *^nvz.*e^ 


-ccsit 

e 


)))<w<Vd' 


■ 
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c(v.  t  ^ 

f  — » 


oO 


i-  ^  \  CASl't  U 


<  ^*>jl=  ^£„£-e 


mKT 


2  dt^'^Cv 


-Lsi'L*  CS/t1  (( 

-,viwo<(ln.l-TwV  "V’^'Xv-X^e  +  ^joi1  U 


+^«'V  -  '^0  e  +  ^  2  (Xl(V.t+€  ~ 


^*n'**£  “  ^VNW^e  +^'A-\V€:  ^  “^C^-v\_ 


z-'-e 


^^-vxtz-e^  +  ^Vj-6  ■  ~  Tn-v-e  **  *VV\-€ 

+  '  Z^W  *  ^tuz-  +T1(V_JS)  e  -  (  1^  *  Z-Ivx 


+  +"^*1+2.^^  “  ^CC  2.  CXy^-Xy^.^  —  £h.t  -V  -^-VWX^)  *— 


USlT1 


- <  *C  tV4  -x^  -x** + Y)  4  ,yl vs  **  ( (■  S*- 


-xj) 


CSl't 

e 


I 


+  (  &  Ux*.z  -XA-X^V^-  '•Xa. 


4  CXv^^+X^_^)  -  aCt^  +  I^i  —  'Z.CX^^* 


-IS/t' 
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+  (  ZX-fe  'X-w  +  -U-.-£  +!»«.., 


+  +  Z"^V2-e  "  ^W-a-fe  e- 


t€Sl't 


+  e  V 


~  tYvci*  ^  2”^”n+e  ”  ^~*V"'+€  '*-^Wt4N)  —  ^X^^  -vTl 


'ma-vg- 


+  2 


twi+fc  *  Xywz+€^)  ~XIn+J^fe  “  ^vv-^+t')) 


cC^-^^Vv  ^  <14  vlv\a^H 

<  io  »  ^  E'xE^t  e.  S  ) )  e.  f«.fe*‘y|,1(tIM 


mvcT 


.  .  -taet1  .  Is^t1  f  __  __  _  v  - 

-Tx„2  -Xy^  6-  *  L^f^\Vxe'  '  ^*^Avi  ‘^“^2  1  *  c^ot 


*  I 

((  ^(X^-X*.^  +  Xy^-Xv^x  A'Xyy  e  +-  Cx^-'ZX^ 


~Xyv_!t  A'Xy^  +  ^^A-S  C.  At  C  y_  ,-2  C  X>4  .  “  ^Wu-Vt^ 


VTiol*1  V  ^-14€ 


_  Xy\-l4e  '  Xn+l+E^  +  ^  ^  -^VV-^e'X-^Z-Vfe^)  '’Xvv-i*.^  ' 


tm1 


+  ^  C  Xn^-X^  -  sCX^-X^  -  4CXa_^-XVV4^ 

*  x. .  .  -T.„  o  ?***'  )  + 16^^  (  Uc  xn.v 


'A-j-e  *+i-€ 


~Xa_t  A-X^^e  A'C  zCX^-Tv^.^  -X^^ArX^^ 


-ia^‘ 

e 


—  £  (  C  C) T^  +■  aCXk-^-^Xvi^^  -T*^-  X^_x^  e* 


iesix1 


•  -f 


. 


-27- 


-i€>Sl-t! 

e 


i 


+  <*  ^^^15  (  (■  ~  ~  ^Xyv-v'^'  -^W-x  ~  “X^  ^X*_x 


mjt* 


‘uSlt1 


+  l^(XA„r  -  Xn+l^  -  £3^  *  +  ^Xy^  +  VUX^  -v  i Xv\-4^)^ 


+  C  “  ^  Xn_v  -  ^^n-2.  "Xv\  '  ^Xvwi  +  '^-Wi  -  V^X\v^ 


"■  ^X^  ^  ^*^VW\  '*’’  ^  ^  '^'Vv+'Z.  "*  6  -L 


vws 


+  e 


-ujv* 


c 


^  G  ^  md1^ n-\-t '”^-n+i-ei  "  5  (x^_v_t  -T  n+\-^l  ^  -X^.^ 


te&V 


X a-x * t  ^X^.^e.  +  t  (  mdC1  ^  Xv-wt  vstxn_\^ 

"  X  »wi-vfc ^  ^  ^XvUfe  - Xfl.y^g'')  ^X^^ +t  e.  ^  ^  ^C>i^  <k*X  A/t* 

The  integration  with  respect  to  ^  can  now  be  performed 
in  equations  2.42-2.44.  Two  types  of  integrals  involving  ^ 
are  encountered  in  equations  2.42-2.44. 


1.  Integrals  of  the  Form 


fe>6 

A*(  H 

O 

where  (  )  is  independent  of  ^  If  is  assumed  to  have 

a  small  positive  imaginary  component  then  this  integral  can 


cC\vQ.*vx\  -  K^O't 


be  evaluated  to  give 


v\a  +wt  - 


2.  Integrals  of  the  Form 

«<> 


krt  C  )'te. 


Integration  by  parts  give  this  integral  the  value 

-  1 _ 

(VUI.  tUi,+  \<VlrV^ 

The  results  in  1  and  2  can  be  expanded  to  give 

it  ^ 


^>31.  tVO  4  V 


and 


-  1 


**  C.  ^  I  ^  r 

V.  WJl+v&t  J 

V  |  l  4-  r 

3.VVC>^  l  J 


(2.45) 


(2.46) 


^vvst-t  -v  (^rt.s.v'C>\S" 

This  expansion  is  valid  providing  is  always  sufficiently  removed 

from  the  cyclotron  frequencies  of  the  ions  of  electrons.  The  approxi¬ 
mation  ^ ^  \ Vx)^\  is  equivalent  to  the  assumption  that  the 

thermal  velocity  over  the  doopler-shif ted  phase  velocity  is  much  less 
than  1.  By  integrating  2.42-2.44,  rearranging  the  terms  in  summations 
and  finally  collecting  terms  of  common  power  in  . 


«>o 
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•2Ljl 


Yt\K~V 


2  i. 
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To  complete  the  calculations  of  the  second-order  currents 


induces  in  the  plasma  by  high  frequency  incident  fields, 
the  integration  over  and  must  be  performed.  All  of 
the  integrals  remaining  to  be  evaluated  are  of  the  form 

“  -  vy>S-^- 


ff  P  cUi+tvs 
t  =  ^  )) Vl  e  **■' 


*ZKV 


"fir  "c*  ° 


where 


P=  0, 1,1,3 

In  order  to  demonstrate  the  convergence  of  Fp  the  square 
in  the  exponent  may  be  completed  and  the  order  of  integration 

•j 

inverted  to  yield  -^KT*'2  ?  -  m  Ns  +  ') 


F  =  Vl 
f  ftf  o 


s 


Art'  e 


2 .xn 


[ 
—  04 


e  V/  AN- 


In  this  form  it  is  easy  to  see  that  Fp  is  finite  and  converges 
rapidly  for  large  values  of  and  •  The  integration 

of  2.47-2.49  can  be  carried  out  to  yield 


(2.50) 


<  *  p  e"  ^  4  <A*'p>  ¥ 

T1‘  1  V\--«4 


(2.51) 


p  s'  v 


(2.52) 
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In  the  regime  under  consideration  in  this  thesis  Xv>l 


Bv  using  this  assumption,  the  Bessel  Function  of  imaginary 
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argument  can  be  expanded  as 

Tn(\V*e  /  \[ 

Through  the  use  of  this  expansion 

I>l  ~  -  #3,  =  ^iz  O 


Avcre 


VY\  ©lZV*TC\ 


=  -  2VC.T  e. 


Fn  cA.2-  x/arTX 


In  this  order  of  approximation 
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2.5  Evaluation  of 


L 


is  given  from  2.18  by 


* 


-  Zee 

3z  \r\  ~oA 


\  Al' 

L  c  av' 


(2.56) 


To  evaluate  ^  >  i  must  first  be  calculated.  Replacing 


^  T^> 

the  subscript  1  by  2  in  equation  2.10,  and  taking  the  complex 


conjugate  yields 


i: 


~  -  zee 


{(C*  li"  «Lt" 

-4  ^ 


(2.57) 


Using  Maxell's  equation, 


* »  WiB. 
C. 


l<„  * 


alone.  Changing 

from  Lagrangian  to  Eulerian  co-ordinates  then  yeilds 
*.  -  vO2k0  P*4 

i  =  \  e  £0(V>(.<  Co-iSi-X  -  £N '  51x5  At 

Z  rn  *^o 

Using  Maxwell's  equation  to  eliminate  in  terms  of  in 

and  expressing  double  primed  quantities  in  terms  of  primed  quan¬ 


tities  yields 
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£rx.l»3\  'S 


$  ^a'-v 


It  is  now  possible  to  carry  out  the  integration  with  respect 


to  ^  .  The  integrals  involved  are  of  the  form 

„->  c(VtV|.'  -vx^.sO't 

Tz-  l  el  '(At 

T4=  \  e  C  '(A't. 

o 

r”4  WivJLlt 
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O 


where  (.  ^  is  independent  of  *£  .  If  ^  is  assumed  to  have 

a  small  negative  imaginary  component,  then  the  integrals 
can  be  evaluated  to  give 


7. 
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We  can  assymptotically  expand  the  integral  to  get 


T,=  -i-C  vC i  V  V^'  ) 

XT=-U  ■)(  U'<oi^ 

X,=  -(.  ")  ( i  ^ 

Xi-  -(.  'iC'  +  ^'O 

(2.59) 

(^-vsiy-  (w^^sr) 

Using  the  results  in  2.59  and  unravelling  the  Lagrangian  co¬ 
ordinates  in  2.58,  and  finally  collecting  terms  of  the  various 
powers  of  V-*.  ^  ^  ^ecomes 

L(.V,  £-wt>  r°  i.CcViVx  +  V>.^+Oi-t')  t  wot* 
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where 
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2.6  Velocity  Moments  of  *n 

•  3t- 

The  results  from  equations  2.34  and  2.41  may  be  used  to 

express  the  velocity  moments  of  v  over  Vx.  as  follows 

T32.  }  J 
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can  again  be  expanded  by  using 

exc^<__ 


Using  equations  2.41  and  rearranging  the  infinite  sums  of 


Bessel  functions,  we  can  write  2.61-2.63  as 
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( Xi-X^  -  ocH  aC  Vs  -XX^X*-XXyO 

(2.65) 
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To  complete  the  calculation  of  the  moments  over  -V  , 

■Viz 

in  the  regime  of  interest  in  this  study,  the  integations  over 
and  'Vi*  are  performed  and  TnC.V'S  is  replaced  by  its  large 
argument  asymptotic  expansion(  XAA  ->  .  The  integrals 

over  and  rt‘  are  of  the  standard  formf^Fj  . 


(2.67) 
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2.7  Evaluation  of  Driving  Currents 

Equations  2.53-2.55  and  equations  2.67-2.69  can  be  now 
used  to  express  the  total  second-order  driving  current  as 

?;A  =  £ £*>  -  < *  (2-70> 
where  k*  1,2.  . 

Equation  2.73  can  be  expressed  as 

3*4  =  (2.71) 

where  C:r  ref  ers  to  the  components  of  the  current. 
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and 
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2.8  Velocity  Moments  of  J) 

*33 

From  the  equation  2.5,^  is  found  to  be  given  by 
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(2.75) 


The  integration  of  this  equation  is  done  in  Stix,  and  the 
results  that  are  applicable  to  this  case  are 
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e  S  fVvv\ 
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J^Note;  the  results  given  by  equation  2.76-2.78  are  obtained 
from  Stix  with  a  Maxwellian  velocity  distribution  assumed 
in  the  z-direction  and  approximating^^^  by  * ts  asymptotic 


value 
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Following  the  analysis  in  Stix,  the  expectation  value  of  the 
velocity  over  the  second-order  induced  distribution  function 
can  be  expressed  as 


<*,$*>*  A-  tfc*  ( 

B0 

where  the  various  components  of  the  mobility  tensor  are  given 
by 
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(2.80) 


(2.81) 


All  other  components  of  the  mobilitv  tensor  are  zero. 

(In  the  limit  ^  ) 

\TzrtX 

2.9  Calculation  of  Total  Current 

Using  2.79,  the  second-order  induced  current  can  be  ex 
pressed  as 


The  total  second-order  current  for  each  species  is  then 


given  by 

and  the  total  current  is  given  by 
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2.10  Solution  of  Second-Order  Fields 

Once  the  second-order  currents  have  been  calculated. 

Maxwell’s  equation  can  be  solved  for  the  second-order  field 
quantities . 

T?  x£.  =  -  et  =  J-3%*  *  ATI  3* 

"  C  5^  "  L  (2.85) 


Fourier  analysing  2.85  and  using  the  refractive  index 
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In  tensorial  form,  this  equation  becomes 
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By  using  2.90  in  2.89,  the  electric  fields  may  be  obtained 
in  terms  of  the  driving  fields. 


2.11  Calculation  of  Energy  Absorption 

The  power  absorbed  by  the  second-order  induced  currents 
can  be  calculated  by  using  equation  2.89  for  the  second-order 
electric  field  and  by  using  equation  2.83  for  second-order 
induced  currents.  The  power  absorbed  by  the  ions  is 


ATKSUr 

The  power  absorbed  by  the  electrons  is  given  by 


(2.91) 


-  ecVi'ovv, 
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and  the  total  power  absorbed  is 
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CHAPTER  3  DISCUSSION  OF  RESULTS 


This  chapter  includes,  as  part  of  a  discussion  of  the 
feasibility  of  ion  cyclotron  heating: 

1.  A  comparison  of  collisional  versus  collisionless  energy 
absorption  of  the  difference  frequency  wave. 

2.  A  comparison  between  first  and  second-order  energy  absorption 
by  the  ions. 

3.  A  complete  analysis  of  the  dependence  of  absorbed  power 
on  number  density,  temperature  and  fluctuations  in  magnetic 
field,  incident  angle,  and  difference  frequency. 

4.  Discussion  of  results  and  suggestions  for  further  research. 


3.1  Comparison  Between  Collisionless  and  Collisional  Absorption 

39 

The  Stix  formula  for  modifying  the  (dimensionless) 
mobility  tensors  to  include  collional  effects  are 


s*e 

54  ”  (3.i) 

A,  ~T 

The  correction  to  the  total  power  absorbed  is  then  found  by 
adding  the  corrections  to  the  mobility  tensors  and  using  this 


new  mobility  tensor  to  evaluate  the  power  absorbed  by  the  ions. 
Calculations  of  the  correction  to  the  power  absorbed 
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bv  collisional  effects  were  carried  out  using 
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The  correction  to  the  power  absorbed  was  found  to  be  less  than 
8%.  Varying  the  nominal  parameters  (i.e.  temperature,  magnetic 
field,  mixing  angle  and  difference  frequency)  placed  and  upper 
bound  of  10%  to  the  ratio  of  collisional  to  collisionless  ion 
power  absorption.  Due  to  their  relative  unimportance,  in  the 
temperature  and  density  range  considered  in  this  study  (for  temp- 
erature  between  £x.io  Kto  b  *  \o  Kand  number  densities  between 
K)U/c.c  .and  )  collisional  corrections  will  be  ignored  in 

the  remainder  of  this  chapter. 

3.2  Comparison  of  First  and  Second  Order  Energy  Absorption 
From  the  Incident  Field 

A  cold  plasma  analysis  gives  the  power  absorbed  by  the  ion  and 
electrons  from  the  incident  wave  through  collisional  damping  as 

VJi*  -  See  I  '^\y\ 

For  the  case  under  consideration  vS2.„>  .Sethis 

reduces  to 

.  S,e  'E,\a 


(3.2) 

} 

expression 

(3.3) 
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The  collision  freqency  used  in  3.2  can  be  estimated  using  the 
Spitzer  formula 

/ (jCT^  ( 3  • 4) 

Computations  to  compare  the  first-order  and  second-order  energy 
absorption  were  carried  out  using  an  I.B.M.-360  for  the  following 
nominal  values. 
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The  resulting  first-order  powers,  using  equations  3. 3-3. 5,  are 

U^e  ^  I  x  E^2* 


Ui1***  5.AxlO^( 


“A  -2. 


The  total  power  absorbed  is 


*m  U  E 


\ 


The  power  absorbed  from  the  difference  frequency  harmonic 
[equations  2.91  and  2.92]  for  the  parameters  given  in  3.5  is 
given  by 


=  2*1  x\o  x  U, 


-\A  ^  4 


and 


e\ec-Vv  ot\ 


(3.6) 


(3.7) 


(3.8) 


uaxurl.6,4 

Comparison  of  first  and  second  order  ion  power  absorption  yields 
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Therefore,  for 


by  the  ions  is  larger  than  the  first-order  power  absorbed  by 
the  ions.  This  means  that  the  non-linear  heating  takes  over  from 


linear  heating  when  fc.  > 


3.3  Dependence  of  Absorbed  Power  on  Fluctuations  in  Plasma 
and  Experimental  Parameters 

Equations  2.91  and  2.92  were  used  to  calculate  the  depen¬ 
dence  of  the  absorbed  power  on  plasma  and  experimental  para¬ 
meters.  The  graphs  shown  in  figures  1-4,  give  the  results  of 
computations  carried  out  using  an  I.B.M.-360  and  the  following 
parameters 


Yi  ~  lo  /  ortr 

r  *  'o'*  °  vc. 

6=  10- * 


In  figure  5,  the  nominal  values  used  were  the  same  as  above 


except  for  Q-  ,  which  is  set  equal  to  20°. 


•M 
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The  nominal  value  of  er  in  the  computations  was  chosen  to  allow 
a  larger  range  of  temperature  variation.  The  magnetic  field 
and  mixing  angle  were  chosen  to  ensure  the  convergence  of 
to 


Magnetic  Field  in  Gauss 
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Ion  Power  =  (  )*  lo'W*  /em'/sec. 

Total  Power  =  C  \o"'\  E,4e.v^s/c.rWsec. 
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Power  Versus  Magnetic  Field  Fluctuation 


-56- 


Ion  Power  = 
Total  Power 


xvo' x  fct^*|c rv\?/s€C. 
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Power  Versus  Frequency  Fluctuations 
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Number  Density  =  Particles/cm. 
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Power  Versus  Number  Density 
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Fig. 


Temperature 
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Ion  Power  =  <  x 

Total  Power  =  C  U  Vo""4*  €^«v^|cw^  |««c. 
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The  computations  oerformed  show  that,  as  \  gets  smaller, 
the  power  absorbed  from  the  second-order  fields  increases. 

This  fact  is  born  out  by  observing  figures  1,  4,  and  5.  QNote; 

A  -  Kx,  KT  J  (  Xn  figure  1,  we  note  the  increase  in  absorbed 
power  with  increasing  magnetic  field.  In  figure  5,  we  note  a 
large  increase  in  the  absorbed  power  as  the  temperature  goes 
down  and  in  figure  4,  we  note  an  increase  in  absorbed  power 
with  a  decrease  in  the  mixing  angle,  -(i.e.  decreasing  Kx  )  • 
Decreasing  corresponds  to  increasing  the  perpendicular 

wavelength  of  the  mixed  wave.  The  large  perpendicular  wave¬ 
length  can  be  achieved  in  two  ways: 

1.  Decreasing  the  mixing  angle 

2.  Decreasing  the  frequency 

For  this  type  of  wave  mixing  phenomena,  the  absorbed  power 

C  4 

goes  as  .  This  would  indicate  that  the  best  way  to  in- 

crease  the  perpendicular  wavelength  would  be  to  decrease  the 
frequency.  (Keeping  in  mind  that  must  be  larger  than 

to  ensure  good  penetration  of  the  incident  fields  into  the 
plasma) . 

The  graphs  plotted  in  figures  1-5  show  the  heating  scheme 
is  relatively  insensitive  to  fluctuations  in  nominal  parameters. 
The  relative  insensitivity  of  the  absorbed  power  to  fluctuations 
brings  the  scheme  under  study  with  the  realm  of  feasibility  as 
a  plasma  heating  method.  For  electric  fields  of  lA\toS 

statvolts/cm.  difference  frequency  heating  takes  over  from  col- 
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lisional  first-order  heating.  Since  electric  fields  of  X\0 

statvolts/cm.^^have  been  achieved  in  high  intensity  lasers,  it 
is  now  possible  to  conduct  an  experiment  to  detect  heating  due 
to  the  difference  frequency  harmonic. 

The  graphs  in  figures  1-5,  show  that  the  power  absorbed 
from  the  difference  frequency  wave  goes  mainly  to  the  electrons. 
This  observation  leads  to  the  conclusion  that  an  analysis 
should  be  carried  out  in  the  regime  where  the  difference  fre¬ 
quency  is  tuned  to  the  electron  cyclotron  frequency.  This  tuning 
to  the  electron  cyclotron  frequency  should  enhance  the  absorbed 

power  by  a  factor  at  least  as  large  as  fill  .  From  the  graphs 

me. 

in  figures  1-5,  we  notice  that  the  absorbed  power  increases  with 
decreasing  \  .  This  observation  leads  to  the  conclusion  that 

the  analysis  done  in  this  thesis  should  be  modified  such  that  the 
asymptotic  expansion  of  the  Bessel  Function  for  small  arguments 
could  be  used  instead  of  the  larger  argument  expansion  used  in 
this  thesis. 

Consideration  could  also  be  given  to  the  mechanism  by 
which  the  difference  frequency  wave-packet  created  in  a  semi¬ 
infinite  plasma  diffuses  in  the  plasma.  This  study  could  have 
practical  implications  since  in  order  to  achieve  fusion  (fusion 
here  refers  to  the  generation  of  useful  energy  from  a  fusion 
reaction)  large  volumes  (large  compared  to  the  spot  size  of  a 

4 

laser  beam)  of  plasma  would  have  to  be  heated  to  ensure  that  a 
fusion  reaction  could  be  maintained.  Since  the  heating  should  be 
non-local,  the  diffusion  time  of  the  wave-packet  can  be  an  impor- 
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tan  t  parameter  in  the  use  of  this  mechanism  for  plasma  heating. 

Note;  fusion  reactors  where  the  volume  of  the  fuel  (usually  in 
the  form  of  a  solid  pellet)  is  of  the  order  of  the  spot  size  of 
a  laser  have  real  possibilities  but  are  not  included  in  this 
discussion. 

The  effects  of  difference  frequency  heating  on  velocity- 
space  instabilities  could  also  be  considered.  i.e.  Does  the 
difference  wave  give  rise  to  large  velocity-space  anisotropies, 
and  if  so  how  do  these  anisotropies  affect  the  macroscopic 
stability  of  the  plasma?  A  more  complete  theory  (for  law  values 
of  T  collisional  heating  becomes  more  important)  would  also 
include  a  fuller  consideration  of  collisional  damping  of  the 
difference  frequency  wave. 
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